Exact Multifractal Exponents for Two-Dimensional Percolation 
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Abstract. The harmonic measure (or diffusion field) near 
a critical percolation cluster in two dimensions (2D) is con- 
sidered. Its moments, summed over the accessible external 
hull, exhibit a multifractal (Mf) spectrum, which I calculate 
exactly. The generalized dimensions D (n) as well as the Mf 
function / (a) are derived from generalized conformal invari- 
ance, and are shown to be identical to those of the harmonic 
measure on 2D random or self-avoiding walks. An exact ap- 
plication to the impedance of a rough percolative electrode 
is given. The numerical checks are excellent. Another set of 
multifractal exponents is obtained exactly for n independent 
self-avoiding walks anchored at the random fractal boundary 
of a percolation cluster. 

PACS numbers: 05.40.Fb, 05.45. Df, 64.60.Ak, 41.20.Cv 



Percolation theory, whose tenuous fractal structures, 
called incipient clusters, present fascinating properties, 
has served as an archetypal model for critical phenomena 
0. The subject has recently enjoyed renewed interest: 
the scaling (continuum) limit has fundamental proper- 
ties, e.g., conformal invariance, which present a mathe- 
matical challenge [||,|],Q. Almost uncharted territory in 
exact fractal studies is the harmonic measure, i.e., the 
diffusion or electrostatic field near an equipotential ran- 
dom fractal boundary, whose self-similarity is reflected 
in a multifractal (Mf) behavior of the harmonic measure 

Mf exponents for the harmonic measure of fractals are 
especially important in two contexts: diffusion- limited 
aggregation (DLA) and the double layer impedance at 
a surface. In DLA, the harmonic measure actually de- 
termines the growth process and its scaling properties 
are intimately related to those of the cluster itself Q. 
The double layer impedance at a rough surface between a 
good conductor and an ionic medium presents an anoma- 
lous frequency dependence, which has been observed by 
electrochemists for decades. It was recently proposed 
that this is at heart a multifractal phenomenon, directly 
linked with the harmonic measure of the rough electrode 
. In both the above contexts, percolation clusters have 
been studied numerically as generic models. 

In this Letter, I consider incipient percolation clusters 
in two dimensions (2D), and determine analytically the 
exact multifractal exponents of their harmonic measure. 
I use recent advances in conformal invariance (linked to 
quantum gravity), which allow for the mathematical de- 
scription of random walks interacting with other random 



fractal structures, such as random walks an( ^ se ^~ 

avoiding walks |L0| . A further difficulty here is the pres- 
ence of a subtle geometrical structure in the percolation 
cluster hull, recently elucidated by Aizenman et al. . 
Excellent agreement with decade-old numerical data is 
obtained, thereby confirming the relevance of conformal 
invariance to multifractality; the exact prediction for the 
anomalous exponent of a percolative electrode given here 
also corroborates the multifractal nature of the latter. 
As an illustration of the flexibility of the method, I also 
give the set of exact multifractal exponents correspond- 
ing to the average nth moment of the probability for a 
self-avoiding walk to escape from a percolation cluster 
boundary. 

Consider a two-dimensional very large incipient cluster 
C, at the percolation threshold p c . Define H (w) as the 
probability that a random walker (RW) launched from 
infinity, first hits the outer (accessible) percolation hull 
Ti(C) at point w e H(C). We are especially interested in 
the moments of H , averaged over all realizations of RW's 
and C 



(1) 



where n can be, a priori, a real number. For very large 
clusters C and hulls TL (C) of average size R, one expects 
these moments to scale as 



Z n » (a/R) T 



(n) 



(2) 



where a is a microscopic cut-off, and where the mul- 
tifractal scaling exponents r (n) encode generalized di- 
mensions D (n), t (n) — (n — 1) D (n) , which vary in a 
non-linear way with n |^2[ , p^| , |l4| , p^[ . Several a priori re- 
sults are known. D(0) is the Hausdorff dimension of the 
support of the measure. By construction, if is a normal- 
ized probability measure, so that r(l) = 0. Makarov's 
theorem here applied to the Holder regular curve 
describing the hull |l7| ], gives the non trivial informa- 
tion dimension r'(l) = D (1) = 1. The multifractal for- 
malism ]i^ , [j~3| , |T4] , ^5t further involves characterizing sub- 
sets Ha of sites of the hull Tt by a Lipschitz-Holdcr 
exponent a, such that their local H-measure scales as 
H (w e H a ) w (a/R) a . The "fractal dimension" / (a) of 
the set H a is given by the symmetric Legendre transform 
of r (n) : 



dr_ 

dn 



(n) , r(n)+f (a) = an, 



da 



a). (3) 
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Because of the ensemble average (|l|), values of / (a) 
can become negative for some domains of a [|l8| . 

This Letter is organized as follows: I first present in 
detail the findings and their potential physical signifi- 
cance and applications, before proceeding with the more 
abstract mathematical derivation. 

My results for the generalized harmonic dimensions for 
percolation are 



D (n) = - 
V ' 2 



V24n + 1 



e [-i+oo), (4) 



valid for all values of moment order n, n — gj. The 



Legendre transform (||) of r (n) = (n — 1) D(n) reads 
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FIG. 1. Universal harmonic multifractal dimensions D(n), 
and spectrum f(a) of a 2D incipient percolation cluster, com- 
pared to numerical results by Meakin et al. (in red). 

Figure 1 shows the exact curve D (n) (|4|) together with 
the numerical results for n € {2,..., 9} by Meakin et 
al. showing fairly good agreement. The slight up- 
wards move from the theoretical curve at high values of 



n suggests a difference between annealed and apparent 
quenched averages, as in the DLA case [ p0[ . 

The first striking observation is that the dimension of 
the support of the measure D (0) ^ Du, where £>h = \ 
is the Hausdorff dimension of the standard hull, i.e., the 
outer boundary of critical percolating clusters In 
fact, -D(O) = | is the dimension Z?ep of the accessible ex- 
ternal perimeter p2] , |il"| , the other hull sites being located 
in deep fjords, which are not probed by the harmonic 
measure. Its exact value Z?ep = | has been recently 
derived in terms of relevant scaling operators describing 
path crossing statistics in percolation pj . In the scal- 
ing continuous regime of percolation, the fjords do close, 
yielding a smoother (self-avoiding) accessible perimeter 
of dimension |. This is in agreement with the instability 
phenomenon observed numerically on a lattice: removing 
the fjords with narrow necks causes a discontinuity of the 
effective dimension of the hull from Dn ~ | to -Dep — I , 
whatever microscopic restriction rules are choosen p2[ . 
In other respects, a 2D polymer at the Q-point is known 
to obey exactly the statistics of a percolation hull p3[ , 
and the Mf results (^j^) therefore apply also to that case. 

An even more striking fact is the complete identity of 
Eqs. (Q|]) to the corresponding results both for random 
walks and self-avoiding walks (SAW's) poj . In particu- 
lar, D (0) = | is the Hausdorff dimension of a SAW, com- 
mon to the external frontier of a percolation hull and of 
a Brownian motion |^||. Seen from outside, these three 
fractal curves are not distinguished by the harmonic mea- 
sure. As we shall see, this fact is linked to the presence 
of a universal underlying conformal field theory with a 
vanishing central charge c = 0. 

The singularity at a = 4 in the multifractal function 
f(a) is due to points on the fractal boundary where the 
latter has the local geometry of a needle. Indeed, by ele- 
mentary conformal covariance, a local wedge of opening 
angle 9 yields an electrostatic potential, i.e., harmonic 
measure, which scales as H(R) ~ R~t> ~ R~ a , thus, 
formally, 9 = — , and 9 — 2tt corresponds to the lowest 
possible value a = ^ . The linear asymptote of the / (a) 
curve for a — > +oo, / (a) ~ — ^ corresponds to the low- 
est part n — ► n* = — ^ of the spectrum of dimensions. 
Its linear shape is quite reminiscent of the case of a 2D 
DLA cluster [M . Define Af (H) as the number of sites 
having a probability H to be hit. Using the Mf formalism 
to change from variable H to a (at fixed value of a/R), 
shows that Af (H) obeys, for H — > 0, a power law behav- 
ior with an exponent t* = 1+ lim —f (a) = 1 + n*. 

a — > + oo a 

Thus we predict 

l^o«#- r ,r* = -. (7) 

This t* = 0.95833... compares very well with the result 
t* = 0.951 ± 0.030, obtained for 10~ 5 < H 10~ 4 

Let us consider for a moment the different, but related, 
problem of the double layer impedance of a rough elec- 



Af(H) 
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trode. In some range of frequencies w, the impedance 
contains an anomalous "constant phase angle" (CPA) 
term (iuj)~^ , where (3 < 1. It was believed that (3 would 
be solely determined by the Hausdorff dimension D (0) 
of the electrode surface. From a natural RW represen- 
tation of the impedance, a different scaling law was re- 
cently proposed: (3 — -^j. (here in 2D), where D (2) 
is the multifractal dimension of the H-measure on the 
rough electrode 0. In the case of a 2D porous percola- 
tive electrode, our results (Q) give D (2) = yi, D (0) = |, 
whence (3 = = 0.6875. This compares very well with 
a numerical RW algorithm result [^5|, which yields an 
effective CPA exponent (3 ~ 0.69, nicely vindicating the 
multifractal description 0. 

Let me now give the main lines of the derivation of ex- 
ponents D (n) by generalized conformal invariance. We 
focus on site percolation on the 2D triangular lattice; by 
universality the results are expected to apply to other 2D 
(e.g., bond) percolation models. The boundary lines of 
the percolation clusters, i.e., of connected sets of occupied 
hexagons, form self- avoiding lines on the dual hexagonal 
lattice. They obey the statistics of loops in the O (N = 1) 
model, where N is the loop fugacity, in the so-called "low- 
temperature phase" , a fact we shall recover below |2^| . 

By the very definition of the H-mcasure, n indepen- 
dent RW's diffusing away from the hull give a geometric 
representation of the n th moment H n , for n integer. The 
values so derived for n S N will be enough, by convexity 
arguments, to obtain the analytic continuation for arbi- 
trary n's. Figure 2 depicts n independent random walks, 
in a bunch, first hitting the external hull of a percolation 
cluster at a site w — (•) . As explained in ref. |0|, such 
a site, to belong to the accessible hull, must remain, in 
the continuous scaling limit, the source of at least three 
non-intersecting crossing paths, noted £3, reaching to a 
(large) distance R. These paths are "monochromatic" : 
one path runs only through occupied (light blue) sites; 
the other two, dual lines, run through empty (white) sites 
jllj . The definition of the standard hull requires only the 
origination, in the scaling limit, of a "bichromatic" pair 
of lines S2 ■ Points lacking additional dual lines are 
not accessible to RW's after the scaling limit is taken, 
because their (white) exit path becomes a strait pinched 
by other parts of the (light blue) occupied cluster. The 
bunch of independent RW's avoids the occupied cluster, 
and defines its own envelope as a set of two boundary 
lines separating it from the occupied part of the lattice, 
thus from S3 (Fig. 2). 

Let us introduce the notation A A B for two sets, A, 
B, of random paths, conditioned to be mutually avoid- 
ing, and A V B for two independent, thus possibly inter- 
secting, sets [ pL0| . Now consider n independent RW's, or 
Brownian paths B in the scaling limit, in a bunch noted 
(yB) n , avoiding a set Sg = (A"P) £ of £ non-intersecting 
crossing paths in the percolation system. Each of the 



latter paths passes only through occupied sites, or only 
through empty (dual) ones. The probability that the 
Brownian and percolation paths altogether traverse the 
annulus T> (a, R) from the inner boundary circle of ra- 
dius a to the outer one at distance R, i.e., are in a "star" 
configuration Si A (V£>)™(Fig. 2), is expected to scale for 
R/a — > 00 as 

V R (S t A n) w (a/R) x{StAn) , (8) 

where we used Si A n = Si A (\/B) n as a short hand no- 
tation, and where x (S£ A n) is a new critical exponent 
depending on £ and n. It is convenient to introduce sim- 
ilar "surface" probabilities Vr (Se A n) ss (a/ R) x ^ Sl!fyn ^ 
for the same star configuration of paths, now crossing 
through the half-annulus T>(a,R) in the half-plane. 
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FIG. 2. An "active" site (•) on the accessible external 
perimeter for site percolation on the triangular lattice. It 
is defined by the existence, in the scaling limit, of £ = 3 
non-intersecting, and "monochromatic" crossing paths £3 
(dotted lines), one on the incipient (light blue) cluster, the 
other two on the dual empty (white) sites. The points are 
entrances of fjords, which close in the scaling limit and won't 
support the harmonic measure. Point (•) is first reached by 
three independent RW's (red, green, blue), contributing to 
H 3 (»). The hull of the incipient cluster (golden line) avoids 
the outer frontier of the RW's (thick blue line). A (random) 
Riemann map of the latter onto the real line R reveals the 
presence of an underlying £ = 3 path-crossing boundary oper- 
ator, i.e, of an L = 4-line, or two-cluster, boundary operator 
in the associated 0(N — 1) model, with dimension in the 
half-plane 2^=3 = x^i^T^L^ = £k=2 ~ 2- Both accessible hull 
and Brownian paths have a frontier dimension |. 

When n —> 0, Vr (Si) ^resp. VR(Si)j is the prob- 
ability of having £ simultaneous monochromatic(non- 
intersecting) path-crossings traversing the annulus in 
the plane (resp. half-plane), with associated exponents 
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X( = x (Se A 0) and xe = x (Se A 0) ||Jll[]. These expo- 
nents have been studied in ref . [[flj , and shown rigorously 
to be actually independent of the coloring of the paths, 
with the restriction in the bulk that there exist at least a 
path on occupied sites and one on dual ones, thus I > 1. 
Here the exponents appear as analytic continuations of 
the harmonic measure ones (ph to n — ► 0, and should 
correspond to the definitions of ref. p| . 

In terms of definition (||), the harmonic measure mo- 
ments (Q) simply scale as Z n » R 2 Vr(Si = 3 A n) [ fTs] , 
which, combined with Eqs. (^ and (p|), leads to 



t (n) = x (S3 A n) 



(9) 



Using the fundamental mapping of the conformal field 
theory (CFT) in the plane R 2 , describing a critical sta- 
tistical geometrical system, to the CFT on a fluctuat- 
ing abstract random Riemann surface, i.e., in presence 
of quantum gravity [p6| , I have recently shown that there 
exist two universal functions U, and V, depending only on 
the central charge c of the CFT, which suffice to generate 
all geometrical exponents involving mutual avoidance of 
random star-shaped sets of paths of the critical system 
0. For c = 0, which corresponds to RW's, SAW's, and 
percolation, these universal functions are: 

U(x) = ±x(l + 2x), V(x) = ±(4x 2 -l). (10) 

with V (x) = U (h (x — ^)). Consider now two arbitrary 
random sets A, 13, involving each a collection of paths 
in a star configuration, with proper scaling crossing ex- 
ponents x (A) ,x(B) , or, in the half-plane, crossing ex- 
ponents x (A) , x (B) . If one fuses the star centers and 
requires A and B to stay mutually avoiding, then the 
new crossing exponents, x (A A B) and x (A A B) , obey 
the star 



10| 



x(AAB) = 2V [If- 1 (x (A)) + U- 1 (x (B))] 
x(AAB) =U [U- 1 (x (A)) + If- 1 (x (B))] , 

where XJ~ X (x) is the inverse function of U 



U- 1 (x) = - (V 24:r + 1 - 1) 



(11) 



(12) 



random paths near their random frontier, i.e., the prod- 
uct of two "boundary operators" on the random surface. 
The latter sum is mapped by the functions U, V, into 
the scaling dimensions in R 2 [ fTpf , The structure thus 
unveiled is so stringent that it immediately yields the 
values of the percolation crossing exponents xe, xe of ref. 
PH ) and our harmonic measure exponents x (Se An) (||). 
First, for a set Se = (AV) e of I crossing paths, we have 
from the recurrent use of <\\A\) 



(13) 



Xe = 2v [eu- 1 ( X1 )] , x t = u [eu- 1 ( Xl )] 



For percolation, two values of half-plane crossing expo- 
nents ie are known by elementary means: ii = l, Xs = 
2 HO. From @ we thus find U' 1 (xi) = \V' X (x 2 ) = 
^U~^x 3 ) — |, (thus xi — ^ which in turn gives 

Xt = 2v {\z) = i (e 2 - 1) , xt = u (\t) = i{i + i). 

We thus recover the identity, previously rigorously estab- 
lished in ref. fPH , of xe — s^l^ -1 ' , xi — x^ZgZp with the 
L-line exponents of the associated O (N = 1) loop model, 
in the "low-temperature phase" . For L even, these expo- 
nents also govern the existence of k = \L spanning clus- 
ters pi] JO] ], with the identity x c k = xi =2 k = ^ (4A; 2 - l) 
in the bulk pl||, and x^ — xe=2k-i — \k (2k — 1) in the 
half- plane ]2l|j28| , |29| . The non-intersection exponents of 
k' Brownian paths are also given by xe, xg for £ = 2k' [0, 
so we observe a complete equivalence between a Brown- 
ian path and two percolating crossing paths, in both the 
plane and half-plane. 

For the harmonic exponents in (|^), we fuse the two 
objects Se and (VS) n into a new star Se A n (see Fig. 
2), and use (|li"|). We just have seen that the boundary 
^-crossing exponent of Si, xe, obeys C/" 1 (xe) = \l- The 
bunch of n independent Brownian paths have their own 
half-plane crossing exponent x((\/B) n ) = nx(B) = n, 
since the boundary dimension of a single Brownian path 
is trivially x (B) — 1 || . Thus we obtain 

x(Se An) = 2V (^ + C/ _1 (n)) . (14) 

Specifying to the case I = 3 finally gives from (JlTf) (|l2|) 



If, on the contrary, A and B are independent and 
can overlap, then by trivial factorization of probabilities, 
x(A\l B) = x (A)+x (B) , and x(A\J B) = x (A)+x (B) 
. The rules (^l|) , which mix bulk and boundary expo- 
nents, can be understood as simple factorization proper- 
ties on a random Riemann surface, i.e., in quantum grav- 
ity H|l0| , or as recurrence relations in M 2 between confor- 
mal Riemann maps of the successive mutually avoiding 
paths onto the line R ||. On a random surface, U^ 1 (x) 
is the boundary dimension corresponding to the value x 
in R x R + , and the sum of U^ 1 functions in Eq. ( O ) 
represents linearly the juxtaposition A A B of two sets of pendent SAW's, namely |T3] x ((VP)") = nx (V) = 



x (S3 A n) = 2 + - (n - 1) + — (V24n + 1 - 5) , 

from which t (n) (|), and D (n) Eq.(@) follow, QED. 

This formalism immediately allows many generaliza- 
tions. For instance, in place of n random walks, one can 
consider a set of n independent self-avoiding walks V , 
which avoid the cluster fractal boundary, except for their 
common anchoring point. The associated multifractal ex- 
ponents x (Si A ((yV) n ) are given by the same formula 
(|l4|), with the argument n in U^ 1 simply replaced by 
the boundary scaling dimension of the bunch of inde- 
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for I = 3. These exponents govern the universal multi- 
fractal behavior of the nth moments of the probability 
that a self-avoiding walk escapes from the random frac- 
tal boundary of a percolation cluster in two-dimensions. 
I thus find that they are identical to those obtained when 
the random fractal boundary is taken as the frontier of a 
Brownian path or a self-avoiding walk fllCfl . 
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